The good, the bad and the boa: An unexpected new species of a true boa revealed by morphological and molecular evidence.

Snakes of the genus Boa are outstanding elements of the New World biota with a broad sociological influence on pop culture. Historically, several taxa have been recognized in the past 300 years, being mostly described in the early days of binomial nomenclature. As a rule, these taxa were recognized based on a suite of phenotypic characters mainly those from the external morphology. However, there is a huge disagreement with respect to the current taxonomy and available molecular phylogenies. In order to reconcile both lines of evidence, we investigate the phylogenetic reconstruction (using mitochondrial and nuclear genes) of the genus in parallel to the detailed study of some phenotypic systems from a geographically representative sample of the cis-Andean mainland Boa constrictor. We used cyt-b only (744bp) from 73 samples, and cyt-b, ND4, NTF3, and ODC partial sequences (in a total of 2305 bp) from 35 samples, comprising nine currently recognized taxa (species or subspecies), to infer phylogenetic relationships of boas. Topologies recovered along all the analyses and genetic distances obtained allied to a unique combination of morphological traits (colouration, pholidosis, meristic, morphometric, and male genitalia features) allowed us to recognize B. constrictor lato sensu, B. nebulosa, B. occidentalis, B. orophias and a distinct lineage from the eastern coast of Brazil, which we describe here as a new species, diagnosing it from the previously recognized taxa. Finally, we discuss the minimally necessary changes in the taxonomy of Boa constrictor complex; the value of some usually disregarded phenotypic character system; and we highlight the urgency of continuing environmental policy to preserve one of the most impacted Brazilian hotspots, the Atlantic Forest, which represents an ecoregion full of endemism.

Trypanosoma sp. infection in Boa constrictor snakes: morphological, hematological, clinical biochemistry, molecular, and phylogenetic characteristics.

There are few reports of Trypanosoma in snakes, as well as little information about its pathogenicity in these animals. Thus, the present study aimed to characterize Trypanosoma found in Boa constrictor snakes, to verify the influence of the parasitism on hematological and clinical biochemistry parameters, and to perform a phylogenetic study of the isolates. Blood samples from sixty-one boas were analyzed for the presence of trypanosomatids and by hematological and clinical biochemistry assays. The flagellates that were found in this analysis were used for cell culture, morphometry, and molecular analysis. Later, molecular typing phylogenetic studies were performed. Nine positive animals (14.75%) were identified by microscopy analysis. The hematological results showed that parasitized animals presented significantly lower levels of packed cell volume, hemoglobin, mean corpuscular volume, and mean corpuscular hemoglobin. In the leukogram, eosinophils and heterophils counts were higher in parasitized animals. Considering the molecular analyses, the isolates presented a higher identity of the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and the 18S small subunit ribosomal RNA (SSU rRNA) gene fragments with Trypanosoma serpentis. The phylogenetic tree, using the GAPDH, clustered all isolates with T. serpentis and Trypanosoma cascavelli. This is the first description of T. serpentis parasitizing boas and of the clinical changes caused by trypanosomatid infection in snakes.

A Multiplex RT-PCR Method for the Detection of Reptarenavirus Infection.

Reptarenaviruses cause Boid Inclusion Body Disease (BIBD), a fatal disease of boid snakes with an economic and ecological impact, as it affects both captive and wild constrictor snakes. The clinical picture of BIBD is highly variable but often only limited. Intracytoplasmic inclusion bodies (IB), which develop in most cell types including blood cells, are the pathognomonic hallmark of BIBD; their detection represents the diagnostic gold standard of the disease. However, IBs are not consistently present in clinically healthy reptarenavirus carriers, which can, if undetected, lead to and maintain the spread of the disease within and between snake populations. Sensitive viral detection tools are required for screening and control purposes; however, the genetic diversity of reptarenaviruses hampers the reverse transcription (RT) PCR-based diagnostics. Here, we describe a multiplex RT-PCR approach for the molecular diagnosis of reptarenavirus infection in blood samples. The method allows the detection of a wide range of reptarenaviruses with the detection limit reaching 40 copies per microliter of blood. Using 245 blood samples with a reference RT-PCR result, we show that the technique performs as well as the segment-specific RT-PCRs in our earlier studies. It can identify virus carriers and serve to limit reptarenavirus spreading in captive snake collections.

Demonstration of Parthenogenetic Reproduction in a Pet Ball Python (Python regius) through Analysis of Early-Stage Embryos.

Parthenogenesis is an asexual form of reproduction, normally present in various animal and plant species, in which an embryo is generated from a single gamete. Currently, there are some species for which parthenogenesis is supposed but not confirmed, and the mechanisms that activate it are not well understood. A 10-year-old, wild-caught female ball python (Python regius) laid four eggs without any prior contact with a male. The eggs were not incubated and, after 3 days, were submitted to the University of Parma for analysis due to the suspicion of potential embryo presence. Examination of the egg content revealed residual blood vessels and a small red spot, indicative of an early-stage embryo. DNA was extracted from the three deceased embryos and from the mother's blood, five microsatellites were analyzed to ascertain the origin of the embryos. The captive history data, together with the genetic microsatellite analysis approach, demonstrated the parthenogenetic origin of all three embryos. The embryos were homozygous for each of the maternal microsatellites, suggesting a terminal fusion automixis mode of development.

Reptarenavirus S Segment RNA Levels Correlate with the Presence of Inclusion Bodies and the Number of L Segments in Snakes with Reptarenavirus Infection-Lessons Learned from a Large Breeding Colony.

Reptarenaviruses cause boid inclusion body disease (BIBD), a fatal disease particularly impacting captive boa constrictor collections. The development of cytoplasmic inclusion bodies (IBs) comprising reptarenavirus nucleoprotein (NP) in many cell types of affected snakes is characteristic of BIBD. However, snakes can harbor reptarenaviruses without showing IBs, hence representing carriers and a potential source of transmission. The RNA genome of reptarenaviruses comprises a small (S) and a large (L) segment, and the snakes with BIBD commonly carry a swarm of reptarenavirus segments. To design sensitive and reliable tools for the diagnosis of reptarenavirus infection in snake colonies, we used metatranscriptomics to determine the reptarenavirus segments present in a large boa constrictor breeding colony. The analysis identified one reptarenavirus S segment and three L segments in the colony. The sequence data served to design real-time reverse transcription-PCR (RT-PCR) targeting the found S segment. This allowed us to identify all infected animals and to quantify the S segment RNA levels, which we found to correlate with the presence of IBs. We further found a positive correlation between the number of L segments and the S segment RNA level, which could suggest that L segment excess also contributes to the IB formation. Information on cohousing of the snakes showed a clear association of reptarenavirus infection with cohousing in general and cohousing with infected animals. Information on breeding and offspring confirmed that vertical transmission occurred. Furthermore, our data suggest that some animals might be able to clear the infection or at least exhibit transient or intermittent viremia. IMPORTANCE Boid inclusion body disease (BIBD) is caused by reptarenavirus infection, and while reptarenavirus nucleoprotein is the main component of the inclusion bodies (IBs) characteristic of BIBD, not all reptarenavirus-infected snakes demonstrate IBs in their cells. Identification of infected individuals is critical for controlling the spread of the disease; however, the genetic divergence of reptarenaviruses complicates reverse transcription-PCR (RT-PCR)-based diagnostics. Here, we tested a next-generation-sequencing-based approach to establish a tailored "colony-specific" set of diagnostic tools for the detection of reptarenavirus small (S) and large (L) genome segments. With this approach, we could demonstrate that an S-segment-specific RT-PCR is highly effective in identifying the infected individuals. We further found the S segment RNA level to positively correlate with the presence of IBs and the number of L segments, which could direct future studies to identify the BIBD pathogenetic mechanisms.

Utility of the burmese Python as a model for studying plasticity of extreme physiological systems.

Non-traditional animal models present an opportunity to discover novel biology that has evolved to allow such animals to survive in extreme environments. One striking example is the Burmese python (Python molurus bivittatus), which exhibits extreme physiological adaptation in various metabolic organs after consuming a large meal following long periods of fasting. The response to such a large meal in pythons involves a dramatic surge in metabolic rate, lipid overload in plasma, and massive but reversible organ growth through the course of digestion. Multiple studies have reported the physiological responses in post-prandial pythons, while the specific molecular control of these processes is less well-studied. Investigating the mechanisms that coordinate organ growth and adaptive responses offers the opportunity to gain novel insight that may be able to treat various pathologies in humans. Here, we summarize past research on the post-prandial physiological changes in the Burmese python with a focus on the gastrointestinal tract, heart, and liver. Specifically, we address our recent molecular discoveries in the post-prandial python liver which demonstrate transient adaptations that may reveal new therapeutic targets. Lastly, we explore new biology of the aquaporin 7 gene that is potently upregulated in mammalian cardiac myocytes by circulating factors in post-prandial python plasma.

A community-science approach identifies genetic variants associated with three color morphs in ball pythons (Python regius).

Color morphs in ball pythons (Python regius) provide a unique and largely untapped resource for understanding the genetics of coloration in reptiles. Here we use a community-science approach to investigate the genetics of three color morphs affecting production of the pigment melanin. These morphs-Albino, Lavender Albino, and Ultramel-show a loss of melanin in the skin and eyes, ranging from severe (Albino) to moderate (Lavender Albino) to mild (Ultramel). To identify genetic variants causing each morph, we recruited shed skins of pet ball pythons via social media, extracted DNA from the skins, and searched for putative loss-of-function variants in homologs of genes controlling melanin production in other vertebrates. We report that the Albino morph is associated with missense and non-coding variants in the gene TYR. The Lavender Albino morph is associated with a deletion in the gene OCA2. The Ultramel morph is associated with a missense variant and a putative deletion in the gene TYRP1. Our study is one of the first to identify genetic variants associated with color morphs in ball pythons and shows that pet samples recruited from the community can provide a resource for genetic studies in this species.

Boid Inclusion Body Disease Is Also a Disease of Wild Boa Constrictors.

Reptarenaviruses cause boid inclusion body disease (BIBD), a potentially fatal disease, occurring in captive constrictor snakes boas and pythons worldwide. Classical BIBD, characterized by the formation of pathognomonic cytoplasmic inclusion bodies (IBs), occurs mainly in boas, whereas in pythons, for example, reptarenavirus infection most often manifests as central nervous system signs with limited IB formation. The natural hosts of reptarenaviruses are unknown, although free-ranging/wild constrictor snakes are among the suspects. Here, we report BIBD with reptarenavirus infection in indigenous captive and wild boid snakes in Costa Rica using histology, immunohistology, transmission electron microscopy, and next-generation sequencing (NGS). The snakes studied represented diagnostic postmortem cases of captive and wild-caught snakes since 1989. The results from NGS on archival paraffin blocks confirm that reptarenaviruses were already present in wild boa constrictors in Costa Rica in the 1980s. Continuous sequences that were de novo assembled from the low-quality RNA obtained from paraffin-embedded tissue allowed the identification of a distinct pair of reptarenavirus S and L segments in all studied animals; in most cases, reference assembly could recover almost complete segments. Sampling of three prospective cases in 2018 allowed an examination of fresh blood or tissues and resulted in the identification of additional reptarenavirus segments and hartmanivirus coinfection. Our results show that BIBD is not only a disease of captive snakes but also occurs in indigenous wild constrictor snakes in Costa Rica, suggesting boa constrictors to play a role in natural reptarenavirus circulation. IMPORTANCE The literature describes cases of boid inclusion body disease (BIBD) in captive snakes since the 1970s, and in the 2010s, others and ourselves identified reptarenaviruses as the causative agent. BIBD affects captive snakes globally, but the origin and the natural host of reptarenaviruses remain unknown. In this report, we show BIBD and reptarenavirus infections in two native Costa Rican constrictor snake species, and by studying archival samples, we show that both the viruses and the disease have been present in free-ranging/wild snakes in Costa Rica at least since the 1980s. The diagnosis of BIBD in wild boa constrictors suggests that this species plays a role in the circulation of reptarenaviruses. Additional sample collection and analysis would help to clarify this role further and the possibility of, e.g., vector transmission from an arthropod host.

Reference genome of the rubber boa, Charina bottae (Serpentes: Boidae).

The rubber boa, Charina bottae is a semi-fossorial, cold-temperature adapted snake that ranges across the wetter and cooler ecoregions of the California Floristic Province. The rubber boa is 1 of 2 species in the family Boidae native to California and currently has 2 recognized subspecies, the Northern rubber boa C. bottae bottae and the Southern rubber boa C. bottae umbratica. Recent genomic work on C. bottae indicates that these 2 subspecies are collectively composed of 4 divergent lineages that separated during the late Miocene. Analysis of habitat suitability indicates that C. bottae umbratica montane sky-island populations from southern California will lose the majority of their habit over the next 70 yr, and is listed as Threatened under the California Endangered Species Act. Here, we report a new, chromosome-level assembly of C. bottae bottae as part of the California Conservation Genomics Project (CCGP). Consistent with the reference genome strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology to produce a de novo assembled genome. The assembly comprises 289 scaffolds covering 1,804,944,895 bp, has a contig N50 of 37.3 Mb, a scaffold N50 of 97 Mb, and BUSCO completeness score of 96.3%, and represents the first reference genome for the Boidae snake family. This genome will enable studies of genetic differentiation and connectivity among C. bottae bottae and C. bottae umbratica populations across California and help manage locally endemic lineages as they confront challenges from human-induced climate warming, droughts, and wildfires across California.

Cytogenetic Analysis of the Members of the Snake Genera Cylindrophis, Eryx, Python, and Tropidophis.

The recent discovery of two independently evolved XX/XY sex determination systems in the snake genera Python and Boa sparked a new drive to study the evolution of sex chromosomes in poorly studied lineages of snakes, where female heterogamety was previously assumed. Therefore, we examined seven species from the genera Eryx, Cylindrophis, Python, and Tropidophis by conventional and molecular cytogenetic methods. Despite the fact that these species have similar karyotypes in terms of chromosome number and morphology, we detected variability in the distribution of heterochromatin, telomeric repeats, and rDNA loci. Heterochromatic blocks were mainly detected in the centromeric regions in all species, although accumulations were detected in pericentromeric and telomeric regions in a few macrochromosomes in several of the studied species. All species show the expected topology of telomeric repeats at the edge of all chromosomes, with the exception of Eryx muelleri, where additional accumulations were detected in the centromeres of three pairs of macrochromosomes. The rDNA loci accumulate in one pair of microchromosomes in all Eryx species and in Cylindrophis ruffus, in one macrochromosome pair in Tropidophis melanurus and in two pairs of microchromosomes in Python regius. Sex-specific differences were not detected, suggesting that these species likely have homomorphic, poorly differentiated sex chromosomes.

Animal Models of Exercise From Rodents to Pythons.

Acute and chronic animal models of exercise are commonly used in research. Acute exercise testing is used, often in combination with genetic, pharmacological, or other manipulations, to study the impact of these manipulations on the cardiovascular response to exercise and to detect impairments or improvements in cardiovascular function that may not be evident at rest. Chronic exercise conditioning models are used to study the cardiac phenotypic response to regular exercise training and as a platform for discovery of novel pathways mediating cardiovascular benefits conferred by exercise conditioning that could be exploited therapeutically. The cardiovascular benefits of exercise are well established, and, frequently, molecular manipulations that mimic the pathway changes induced by exercise recapitulate at least some of its benefits. This review discusses approaches for assessing cardiovascular function during an acute exercise challenge in rodents, as well as practical and conceptual considerations in the use of common rodent exercise conditioning models. The case for studying feeding in the Burmese python as a model for exercise-like physiological adaptation is also explored.

Short '1.2× Genome' Infectious Clone Initiates Kolmiovirid Replication in Boa constrictor Cells.

Human hepatitis D virus (HDV) depends on hepatitis B virus co-infection and its glycoproteins for infectious particle formation. HDV was the sole known deltavirus for decades and believed to be a human-only pathogen. However, since 2018, several groups reported finding HDV-like agents from various hosts but without co-infecting hepadnaviruses. In vitro systems enabling helper virus-independent replication are key for studying the newly discovered deltaviruses. Others and we have successfully used constructs containing multimers of the deltavirus genome for the replication of various deltaviruses via transfection in cell culture. Here, we report the establishment of deltavirus infectious clones with 1.2× genome inserts bearing two copies of the genomic and antigenomic ribozymes. We used Swiss snake colony virus 1 as the model to compare the ability of the previously reported "2× genome" and the "1.2× genome" infectious clones to initiate replication in cell culture. Using immunofluorescence, qRT-PCR, immuno- and northern blotting, we found the 2× and 1.2× genome clones to similarly initiate deltavirus replication in vitro and both induced a persistent infection of snake cells. The 1.2× genome constructs enable easier introduction of modifications required for studying deltavirus replication and cellular interactions.

Phylogeography, historical demography and systematics of the world's smallest pythons (Pythonidae, Antaresia).

Advances from empirical studies in phylogeography, systematics and species delimitation highlight the importance of integrative approaches for quantifying taxonomic diversity. Genomic data have greatly improved our ability to discern both systematic diversity and evolutionary history. Here we combine analyses of mitochondrial DNA sequences, thousands of genome-wide SNPs and linear and geometric morphometrics on Antaresia, a clade of four currently recognised dwarf pythons from Australia and New Guinea (Antaresia childreni, A. stimsoni, A. maculosa and A. perthensis). Our integrative analyses of phylogenetics, population structure, species delimitation, historical demography and morphometrics revealed that the true evolutionary diversity is not well reflected in the current appraisal of the diversity of the group. We find that Antaresia childreni and A. stimsoni comprise a widespread network of populations connected by gene flow and without evidence of species-level divergence among them. However, A. maculosa shows considerable genetic structuring which leads us to recognise two subspecies in northeastern Australia and a new species in Torres Strait and New Guinea. These two contrasting cases of over and under estimation of diversity, respectively, illustrate the power of thorough integrative approaches into understanding evolution of biodiversity. Furthermore, our analyses of historical demographic patterns highlight the importance of the Kimberley, Pilbara and Cape York as origins of biodiversity in Australia.

Phylogenomics, biogeography and taxonomic revision of New Guinean pythons (Pythonidae, Leiopython) harvested for international trade.

The large and enigmatic New Guinean pythons in the genus Leiopython are harvested from the wild to supply the international trade in pets. Six species are currently recognized (albertisii, biakensis, fredparkeri, huonensis, meridionalis, montanus) but the taxonomy of this group has been controversial. We combined analysis of 421 nuclear loci and complete mitochondrial genomes with morphological data to construct a detailed phylogeny of this group, understand their biogeographic patterns and establish the systematic diversity of this genus. Our molecular genetic data support two major clades, corresponding to L. albertisii and L. fredparkeri, but offer no support for the other four species. Our morphological data also only support two species. We therefore recognize L. albertisii and L. fredparkeri as valid species and place L. biakensis, L. meridionalis, L. huonensis and L. montanus into synonymy. We found that L. albertisii and L. fredparkeri are sympatric in western New Guinea; an atypical pattern compared to other Papuan species complexes in which the distributions of sister taxa are partitioned to the north and south of the island's central mountain range. For the purpose of conservation management, overestimation of species diversity within Leiopython has resulted in the unnecessary allocation of resources that could have been expended elsewhere. We strongly caution against revising the taxonomy of geographically widespread species groups when little or no molecular genetic data and only small morphological samples are available.

Revisiting the Karyotype Evolution of Neotropical Boid Snakes: A Puzzle Mediated by Chromosomal Fissions.

The Boidae family is an ancient group of snakes widely distributed across the Neotropical region, where several biogeographic events contributed towards shaping their evolution and diversification. Most species of this family have a diploid number composed of 2n = 36; however, among Booidea families, the Boidae stands out by presenting the greatest chromosomal diversity, with 2n ranging between 36 and 44 chromosomes and an undifferentiated XY sex chromosome system. Here, we applied a comparative chromosome analysis using cross-species chromosome paintings in five species representing four Boidae genera, to decipher the evolutionary dynamics of some chromosomes in these Neotropical snakes. Our study included all diploid numbers (2n = 36, 40, and 44) known for this family and our comparative chromosomal mappings point to a strong evolutionary relationship among the genera Boa, Corallus, Eunectes, and Epicrates. The results also allowed us to propose the cytogenomic diversification that had occurred in this family: a process mediated by centric fissions, including fission events of the putative and undifferentiated XY sex chromosome system in the 2n = 44 karyotype, which is critical in solving the puzzle of the karyotype evolution of boid snakes.

Phylogenomics, Biogeography, and Morphometrics Reveal Rapid Phenotypic Evolution in Pythons After Crossing Wallace's Line.

Ecological opportunities can be provided to organisms that cross stringent biogeographic barriers towards environments with new ecological niches. Wallace's and Lyddeker's lines are arguably the most famous biogeographic barriers, separating the Asian and Australo-Papuan biotas. One of the most ecomorphologically diverse groups of reptiles, the pythons, is distributed across these lines, and are remarkably more diverse in phenotype and ecology east of Lydekker's line in Australo-Papua. We used an anchored hybrid enrichment approach, with near complete taxon sampling, to extract mitochondrial genomes and 376 nuclear loci to resolve and date their phylogenetic history. Biogeographic reconstruction demonstrates that they originated in Asia around 38-45 Ma and then invaded Australo-Papua around 23 Ma. Australo-Papuan pythons display a sizeable expansion in morphological space, with shifts towards numerous new adaptive optima in head and body shape, coupled with the evolution of new micro-habitat preferences. We provide an updated taxonomy of pythons and our study also demonstrates how ecological opportunity following colonization of novel environments can promote morphological diversification in a formerly ecomorphologically conservative group. [Adaptive radiation; anchored hybrid enrichment; biogeography; morphometrics; snakes.].

Genetic diversity in the US ex situ populations of the endangered Puerto Rican Boa, Chilabothrus inornatus.

The Puerto Rican Boa (Chilabothrus inornatus) was placed on the US Endangered Species List in 1970. Progress has been made since to clarify the recovery status of this species, though the design of a new recovery plan must include information regarding genetic variation within and among populations of this species. While measures of genetic diversity in wild populations of this species are finally becoming available, relative genetic diversity represented in ex situ populations is unknown, which hampers efforts to develop an ex situ species management plan. Here, we provide an analysis of genetic diversity in US public and private collections (zoos and breeders) using mitochondrial sequence data and five highly polymorphic nuclear microsatellite loci. We analyzed 50 boas from the US ex situ population and determined overall genetic diversity and relatedness among these individuals. We then compared these data to mitochondrial and microsatellite data obtained from 176 individuals from wild populations across the native range of the species. We found little inbreeding and a large amount of retained genetic diversity in the US ex situ population of C. inornatus relative to wild populations. Genetic diversity in the ex situ population is similar to that found in wild populations. Ours is only the second explicit attempt to characterize genetic diversity at the molecular level in ex situ populations of boid snakes. We anticipate that these results will inform current breeding strategies as well as offer additional information that will facilitate the continuation of ex situ conservation breeding or management in boas.

Species delimitation and systematics of the green pythons (Morelia viridis complex) of melanesia and Australia.

Molecular data sets and the increasing use of integrative systematics is revealing cryptic diversity in a range of taxa - particularly in remote and poorly sampled landscapes like the island of New Guinea. Green pythons (Morelia viridis complex) are one of the most conspicuous elements of this island's fauna, with large numbers taken from the wild to supply international demand for exotic pets. We test hypotheses about species boundaries in green pythons from across New Guinea and Australia with mitochondrial genomes, 389 nuclear exons, and comprehensive assessment of morphological variation. Strong genetic structuring of green python populations and species delimitation methods confirm the presence of two species, broadly occurring north and south of New Guinea's central mountains. Our data also support three subspecies within the northern species. Subtle but consistent morphological divergence among the putative taxa is concordant with patterns of molecular divergence. Our extensive sampling identifies several zones of hitherto unknown biogeographical significance on the island of New Guinea. We revise the taxonomy of the group, discuss the relevance of our findings in the context of Papuan biogeography and the implications of our systematic changes for the conservation management of these taxa.

Cytogenetic Analysis Did Not Reveal Differentiated Sex Chromosomes in Ten Species of Boas and Pythons (Reptilia: Serpentes).

Homologous and differentiated ZZ/ZW sex chromosomes (or derived multiple neo-sex chromosomes) were often described in caenophidian snakes, but sex chromosomes were unknown until recently in non-caenophidian snakes. Previous studies revealed that two species of boas (Boa imperator, B. constrictor) and one species of python (Python bivittatus) independently evolved XX/XY sex chromosomes. In addition, heteromorphic ZZ/ZW sex chromosomes were recently revealed in the Madagascar boa (Acrantophis sp. cf. dumerili) and putatively also in the blind snake Myriopholis macrorhyncha. Since the evolution of sex chromosomes in non-caenophidian snakes seems to be more complex than previously thought, we examined ten species of pythons and boas representing the families Boidae, Calabariidae, Candoiidae, Charinidae, Pythonidae, and Sanziniidae by conventional and molecular cytogenetic methods, aiming to reveal their sex chromosomes. Our results show that all examined species do not possess sex-specific differences in their genomes detectable by the applied cytogenetic methods, indicating the presence of poorly differentiated sex chromosomes or even the absence of sex chromosomes. Interestingly, fluorescence in situ hybridization with telomeric repeats revealed extensive distribution of interstitial telomeric repeats in eight species, which are likely a consequence of intra-chromosomal rearrangements.

Genomic Basis of Convergent Island Phenotypes in Boa Constrictors.

Convergent evolution is often documented in organisms inhabiting isolated environments with distinct ecological conditions and similar selective regimes. Several Central America islands harbor dwarf Boa populations that are characterized by distinct differences in growth, mass, and craniofacial morphology, which are linked to the shared arboreal and feast-famine ecology of these island populations. Using high-density RADseq data, we inferred three dwarf island populations with independent origins and demonstrate that selection, along with genetic drift, has produced both divergent and convergent molecular evolution across island populations. Leveraging whole-genome resequencing data for 20 individuals and a newly annotated Boa genome, we identify four genes with evidence of phenotypically relevant protein-coding variation that differentiate island and mainland populations. The known roles of these genes involved in body growth (PTPRS, DMGDH, and ARSB), circulating fat and cholesterol levels (MYLIP), and craniofacial development (DMGDH and ARSB) in mammals link patterns of molecular evolution with the unique phenotypes of these island forms. Our results provide an important genome-wide example for quantifying expectations of selection and convergence in closely related populations. We also find evidence at several genomic loci that selection may be a prominent force of evolutionary change-even for small island populations for which drift is predicted to dominate. Overall, while phenotypically convergent island populations show relatively few loci under strong selection, infrequent patterns of molecular convergence are still apparent and implicate genes with strong connections to convergent phenotypes.